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Abstract 

Using the 2D Jordan-Wigner transformation we reformulate the square-lattice s = i XY (XZ) model 
in terms of noninteracting spinless fermions and examine the ground-state and thermodynamic properties of 
this spin system. We consider the model with two types of anisotropy: the spatial anisotropy interpolating 
between 2D and ID lattices and the anisotropy of the exchange interaction interpolating between isotropic XY 
and Ising interactions. We compare the obtained (approximate) results with exact ones (ID limit, square- 
lattice Ising model) and other approximate ones (linear spin-wave theory and exact diagonalization data for 
finite lattices of up to N = 36 sites supplemented by finite-size scaling). We discuss the ground-state and 
thermodynamic properties in dependence on the spatial and exchange interaction anisotropies. We pay special 
attention to the quantum phase transition driven by the exchange interaction anisotropy as well as to the 
appearance/disappearance of the zero-temperature magnetization in the quasi-lD limit. 
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1 Introduction 

The s = i Heisenberg model on a square lattice is the well-known basic model in the quantum theory of magnetism 
which became especially attractive after the discovery of high-temperature superconductors pi 0]. Much more 
materials which can be viewed as realizations of the 2D s = \ Heisenberg model are known at present time. Many 
analytical and numerical studies analyzing thermodynamics, spin correlations and their dynamics for the 2D s = h 
Heisenberg model were performed during the last fifteen years. One of the (approximate) analytical approaches 
is based on the 2D Jordan-Wigner fermionization. An extension of the ID Jordan- Wigner transformation for 
higher dimensions was suggested in the late 80s (see also Such kind of approximate treatment was 

applied to the square-lattice quantum spin models in several papers ||, @] (for a brief review see §§). Recently, 
the 2D Jordan-Wigner mapping has been used for the calculation of the magnetization curves of several 2D spin 
systems |{| . The reported studies J|, f| |], @, |§] deal with the s — \ isotropic Heisenberg or the s = \ isotropic XY 
models and the effects of anisotropic XY exchange interaction were not considered. Another intriguing problem 
in the theory of magnetism is the influence of the spatial anisotropy on the properties of quantum spin models 
between the limiting cases of a system of noninteracting chains and of a spatially isotropic 2D system. The issue 
of disappearing of long-range order in ID limit although discussed by several authors (see, for example, the papers 
[ full [il] , [l2[ |l3| on the crossover from one to two dimensions for the s = ^ Heisenberg model) remains ambiguous. 

In this paper we consider the s = | anisotropic XY model on a spatially anisotropic square lattice using the 2D 
Jordan-Wigner fermionization approach in order 1) to examine the effects of exchange interaction anisotropy and 
the effects of spatial anisotropy on the ground-state and thermodynamic properties of the spin model and 2) to 
gain an understanding how well this technique works comparing the approximate results with the exact ones (ID 
limit, the Onsager solution for the square-lattice Ising model). To examine the quality of the 2D Jordan-Wigner 
fermionization scheme in the cases for which exact results are not available we perform linear spin-wave theory 
calculations and finite size scaling analysis of exact diagonalization data. The 2D Jordan-Wigner fermionization 
may be used for a study of the dynamic properties of 2D quantum spin systems jl4| and a first test of this approach 
for calculating the static quantities is strongly desirable. 

There has been a great deal of theoretical work on the 2D s = \ XY model f|, 0, || |§, f|| ||, ||, 
^3| , [24|, |2(J. We use the results reported in these papers in our study for comparison in due course. Thus, in 
Ref. jl7| it was shown how the spin-wave theory can be applied to the quantum XY Hamiltonians. Already the 
very early studies on square-lattice s = \ isotropic XY model indicated the long-range order in the ground state. 
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Later the existence of long-range order in the ground state of quantum XY model on a hypercubic lattice in two 
and higher dimensions was proved rigorously [|l8]. The exact diagonalization study for the isotropic XY model 
on finite square lattices of up to 6 x 6 sites supplemented by finite-size scaling was reported in Ref. . The 
high precision quantum Monte Carlo results presented in Ref. |E| concern with both the properties at zero and 
nonzero temperatures. The zero-temperature properties of the square-lattice s = \ anisotropic XY model were 
studied via spin-wave theory and via series expansions around the Ising limit p3| . The papers |m[ [2^] contain 
the zero-temperature results of the coupled cluster method and the correlated basis function method applied on 
the square-lattice s = ^ anisotropic XY model, respectively. The Green's-function approach was used for the 
calculation of the temperature dependences of the susceptibilities and specific heat of the square-lattice s = j 
isotropic XY model in Ref. j26|. It should be remarked that the mentioned studies refer to the spatially isotropic 
lattices. 

This paper is organized as follows. In Section 2 we present the conventional linear spin-wave theory for 
the square-lattice s = \ XY model. This consideration is a straightforward generalization of Ref. |l7j for the 
anisotropic case and is given here for easy references. In Section 3 we present the 2D Jordan- Wigner fermionization 
treatment of the XY and XZ Hamiltonians. The results derived in this Section are parallel to those obtained 
on the basis of the spin-wave theory and resemble strongly some outcomes of the coupled cluster method and 
the correlated basis function method . In Section 4 we present the finite-size scaling analysis of the exact 

diagonalization data. In Section 5 we discuss the effects of anisotropics on the ground-state and thermodynamic 
properties of the square-lattice s = k XY model. We compare the results of different approaches focusing on the 
validity of the approximate treatment based on the 2D Jordan-Wigner fermionization. Finally, we summarize our 
findings in Section 6. Some preliminary results of this study were announced in the conference papers pTj]. 

2 The model. Linear spin— wave theory 

We consider a model consisting of N = N x N y (N x = N y = y/~N — > oo) spins | on a spatially anisotropic square 
lattice governed by the anisotropic XY Hamiltonian 

n x -i Ny-l 
i=0 j=0 

+J_l ((1 +7K,< J+ i + (i -7)4>^- + i)) 

^ JV X -1 Ny-l 

i=0 j=0 



+ J -i ( s i,j s i,j+i + s i,j s t,j+i + 7 ( s t,j s t,j+i + s i,j s i,j+i))) ■ I 2 - 1 ) 

Here J and J± (= RJ) are the exchange interaction strengths between neighbouring sites in a row (horizontal 
direction) and a column (vertical direction) , respectively, and the parameter 7 controls the anisotropy of exchange 
interaction. Since the considered lattice is a bipartite one, the cases of ferromagnetic and antiferromagnetic signs 
of exchange interactions are related to each other by a simple spin rotation. Therefore, in what follows we assume 
without loss of generality the ferromagnetic sign of the exchange interactions J and J±, i.e., J, Jj_ < 0. We can 
recover the ID limit putting either J± = or J = arriving at a system of independent chains running in horizontal 
or vertical direction, respectively. The case 7 = corresponds to the isotropic XY interaction, whereas the case 
7 = 1 (or 7 = —1) corresponds to the Ising interaction. We are interested in the ground-state and thermodynamic 
properties of spin model (|2.l| ). 

We begin with the linear spin-wave theory following the scheme developed in |Ti|] . The conventional spin- 
wave analysis starts from a ferromagnetic state fully polarized in z direction J|^]. We have to deal with the 
XZ model rather than the XY model. The former Hamiltonian arises from the latter one after the unitary 
transformation s x — > — s z , s y — > s x , s z — > — s y , and hence exhibits the same thermodynamics. (Evidently, within 
approximate approaches in which the x, y spin components are treated differently than the z spin components we 
may expect different results for the XY and XZ models.) The performed rotation is also employed within the 
coupled cluster method p4| or correlated basis function method pq] applied to the quantum XY model (but not 



for exact diagonalization computations). Thus, we consider 

N*-lN y -l 



A=EE (J((l-7K,< +1 , 3 + (1 + 7K^- 



i=0 j=0 

+J X ((1 - 7K J <,+i + (1 + 7)<i<i+i)) 

i=0 j=0 ^ 



1 — 7 

( S i,j S i+lJ + S i,j S i+l,j °i,j s i+l,j ^ s i,j s i+hjJ 



1 — 7 4 

+J± 4 ( s i,j s i,j+l + s i,j s i,j+l + S i,j S i,j+l + s i,j s i,i+l) 

+^ (! + 7) - i) W - 5)) ■ (2-2) 

Employing the Holstein-Primakoff transformation we arrive at the following Hamiltonian of the linear spin-wave 
theory 

# = E E ( 3 ( a tj a t+l,j + a i,j a i+h3 + a i,j a i+l,j + »-,,"/•:.,) 
i=0 j=0 ^ 



. (\ 1 , 1 , 

+■7(1+7) I | - ^ - tj",'. i,,"- lj 
1 — 7 



+J x (l + 7) Q - i"/.,"-., - ^a^+i ai 'i+ 1 )) ' ( 2 - 3 ) 



where the operators afj, fljj obey the Bose commutation rules. After performing the Fourier transformation 

a- ■ - 1 V e^ k " i+k "^ai. i a + - 1 V e ~ i( - k ' i+ ^a + (2 4) 

V iV ^ v ?/ kx ,k y V^x^y kx ,k y 

with fe x = §^n x , n x = -?f, + 1, . . . , ^ - 1, fcj, = j^n v , n y = + 1, . . . , ^ - 1 (JV a , N y are even) 

we come to the following quadratic form 

h = \Y1 W 1 - T) J k - 2 (1 + 7) (J+ J±)) K a k + a-kfl^ k ) + (1 - 7) Jk «a^ k + a k a_ k )) 

k 

+^iV(l + 7)(J + Jx), (2.5) 

J k = J COS fc^ + Jx COS fcj, . 

This quadratic form ( [2.5| ) can be diagonalized by the Bogolyubov transformation 

a k = (coshtik) fe k + (sinhitk) fri" k , 
a+ k = (sinhM k ) 6 k + (coshu k ) bl" k , (2.6) 

where the real function u k — it_ k is given by the equation 

tanh(2u k ) = - (1-7) Jk (2 ?) 

V kj 2(1 + 7) (J + J x )- (I-7) J k 



As a result we get 



H = E n * ( 6+6k + 5) + \n(1 +l)(J + Jx), (2.8) 



= V(l + 7) 2 (J + J±f - (1 ~ 7 2 ) (J + JjJ Jk- 



In the special case 7 = 0, J = Jj_ Eq. ( |2.S| ) reproduces the result derived in Ref. |Ti 
Now we can calculate in the usual way the ground-state energy per site 



I f w dk x dk v „ 3, . . . 



and the zero-temperature (z) magnetization per site 

1 f dk 



2 J ( 2 +> 2 yi - tanh 2 (2u k ) 



+ 1 (2.10) 



(keeping in mind that s 2 



a + a). 



Let us note that the approximate treatment which starts from Eq. (2.3) has destroyed the symmetry with 
respect to the change 7 to —7. Moreover, the consideration is valid only for 7 > otherwise there exist such values 



of k at which the r.h.s. of (2.7) for tanh (2itk) exceeds 1. Evidently, while getting ( [2.3] ) we have presumed the 
long-range order in z direction in the spin space and this assumption becomes wrong if 7 < 0. 



3 2D Jordan— Wigner fermions: XY vs. XZ 

An alternative approach to calculate the thermodynamic quantities of the square-lattice s — i XY model consists 
in the 2D Jordan- Wigner transformation reformulating the problem in fermionic language. First we introduce the 
annihilation and creation operators of spinless fermions 

a *.i = E E Imln( ^ + - j)) d ~f,g d f>9 

transforming fl2.l| ) into 

N X -lNy-l 

i=0 j=0 

+7 (d+ J c- i ^- +a ^d+ +l j + d h3 e i(a '^ +a ^d l+1 ^ 

+7 (dt J e- i ^-^ +a ^d+ ]+1 + ( /,,e'"' ' " d, ,.^)) . (3.2) 

We adopt a mean-field treatment of the phase factors which have appeared in the Hamiltonian (3.2) replacing the 
fermion occupation-number operators by their average value \ and using the same gauge as in Ref. (see also 

Refs. [|[ ||). As a result we arrive at 

N x -iN y -i 



2 

i=0 3=0 



+Jl (d+ - d hj d+ J+1 + 7 (df tj df j+1 - di tj di,j+i))) ■ (3.3) 



This is the only approximate step made to derive the thermodynamic quantities of the spin model (2.1) within the 
2D Jordan- Wigner fermionization approach. Further consideration requires no approximations. Next we perform 
the Fourier transformation 

J. . — 1 \ " pi(*»*+fc»J')-y, , ,1+ — 1 \ " p -Hk x i+k v j) i+ (oa\ 

J ~ /WW ^ iJ ~ /aH\T ^ [ ' 



2 



2 



2 



1, fe 



JVj, 
2 



+ 1, 



2 



1 (JVjb, A y are even) 



As a result Eq. (3.3) becomes 



A LB -iC -£> 
-LB -A L> iC 
iC L> -A -LB 



\ -£) -iC \B A j 



di. 



A = (7j^ cos ky , _B = 7Jisinfc y , C = Jsmk x , D = ~/Jcosk x 



(3.5) 



Here the prime denotes that k in the thermodynamic limit varies in the region — tt < k y < it, — 7t + \k y \ < k x < 



7r — \k y \. The 4x4 matrix in (3.5) can be diagonalized. Its eigenvalues A a (k) are given by 



Ai(k) = 
A 2 (k) = 



(Jj_ cos ky + jj cos k x ) 2 + (J sin k x + 7 Ji sin ky) 2 , 



(Jj_ cosky — jJcosk x ) 2 + (Jsinfc^ — jJx sinfe^) 2 , 



Therefore, the Hamiltonian (3.5) then assumes the form 



. , xVk,, 

"~ a=l ^ 

where ?7k,a, a' 01 ~ 1' 2 are Fermi operators. 

It is easy to calculate now the Helmholtz free energy per site 



A 3 (k) = -A x (k), 
A 4 (k) = -A 2 (k). 



2/' 



In 2 cosh 



£A a (k ) 

2/3 7.,, 2tt 2tt V"" V 2 / ' ~ V 2 

The ground-state energy per site which follows from ( |3.8Q in the zero-temperature limit /3 — > 00 reads 

1 /• dk 



eo 



(27T) 



(Ai(k) + A 2 (k)) 



The entropy and specific heat according to Eq. (3.8) are as follows 



dk 'in f 2 cosh 



in 2 cosh 



/3A 2 (k) 



2 7 (2^) 2 v 2 22 2 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



(3.10) 



c _ 1 

k ~ 2 



dk 

(2^ 



/3Ai(k) 
2 



/3A 2 (k) 
2 



cosh 



cosh 



/3A 2 (k) 



(3.11) 



In ID limit formulas ( |3.6| ), ( ]3 .8| ) - (3.11) yield the long-known exact results for the thermodynamic quantities of 
spin-j anisotropic XY chain. Note, that the symmetry with respect to the change 7 to —7 inherent in the original 



spin model is present in the fermionic description (3.7), (3.6). 

We cannot calculate the magnetization m x (or m v ) by averaging sfj (or sf j) because of the invariance of the 



Hamiltonian ([2.1D when all s 



-sf ■ (or all s\ 



-sf •); this symmetry remains in the approximate treatment 



as well. m x (or m v ) can be obtained as a square root of the two-site correlation function of x (or y) spin components 
taken at sites which are at infinitely large distance. However, to our best knowledge such correlation functions 
have not been examined within the 2D Jordan- Wigner fermionization approach up till now. 

In the second part of this Section we consider the XZ rather than the XY Hamiltonian. We start from the 



Hamiltonian (2.2). Making use of the transformation (3.1) and adopting the mean-field treatment of phase factors 



as used already for obtaining (3.3) we get the following Hamiltonian of interacting spinlcss fcrmions on a square 
lattice 



h = E E 

i=0 j=0 



( d tj d t+i,j + d t,jdi+i,j - <l,. ,'i;. 1 , - <!,.,<!, . 1. 



j d i+i,j - 2 



-Jx 



1-7 



( d tj d t] + l + d tj d i,] + l - d iJ d tj + l _ d i,j d i,j+l.) 

+Jx (i + 7) (<//...</,., - 5) ('/;,., - 5)) ■ 



(3.12) 



We need further approximations because of the interaction terms in Eq. (3.12). We may assume that 



1 



1 



df .dijdr^ .di+ij -> mdf .dij + -'/</,.., + md+d l+hj + -dJ +l j d i+X j - - - 



2 -,., 
1 
2 



d tjdi, J dt :j+1 d l j +1 -> md^dij + ^df^ + md+ j+1 d id+1 + ■zd+ j+1 d i , j + 



1 

2' 



4 
1 

1 — — — m — to, 



(3.13) 



i.e., (ferromagnetic) long-range order (the existence of which in two dimensions is proved rigorously for 7 = at 
zero temperature or for 7 = 1 at nonzero temperature) is imposed. The order parameter is m = m z and its nonzero 
value may be expected for positive 7. For negative 7 one may apply the unitary transformation s x — > s z , s y — > ,s y , 
s z — » — s x coming back to the model with positive 7. Thus the thermodynamics of the model with positive 7 and 
to = m z is the same as of the model with negative 7 and to = m x . These symmetry arguments permit us to find 



the thermodynamic properties of the model (2.2) also for 7 < basing on the results calculated for 7 > 0. 
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Performing further the Fourier transformation (3A) one gets 



( A + M 
-iB 
iC 



2 

iB -iC -X> 

-.4 - X V iC 

V -A + M -iB 

-iC iB A-M J 



^*x i 

^fcx±7T,/Cy ±7T 

\ ^-fc»±7r,-fe s ±7r / 



JV(1 + 7 )(J + J±)m 2 



with 



A = — r— ^- Jj_ cos fcj, , B = — ^-L sm j C = — -^p- J sin k x , V = — — — J cos k x , 

M = 2(1 + 7 ) (J +J ± )m. 



The 4x4 matrix in ( 3.14 ) can be diagonalized. Its eigenvalues A a (k) are given by 



A lj2 (k) = \/A 2 + B 2 + C 2 + V 2 + M 2 ±2^J(AV + BCy + (A 2 +C 2 )M 2 , 

A 3)4 (k) = -A 1)2 (k). 



Therefore, the resulting fermionic Hamiltonian assumes the form 

2 

' A„.fkl 



' / l\ 

h = E A «( k ) ( - o ) - N 0- + 7) ( J + ™ : 



k a=l 



with A Q (k) ( J3.16D - 

The Hclmholtz free energy per site is given by 



2/3 7 (2tt)' 



< lk ^i n f2cosh^iM|+ln[2cosh 



/3A 2 (k) 



+ (l+ 7 )(|J| + |^|)m 2 



and the magnetization m is determined by minimizing (3.18|), i.e., from the following equation 



2(1 + 7) (\J\ + mm =U^ tanh ^ + ^ tanh ^ 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



4 J (2tt) 2 V 9m """" 2 ' 9m """" 2 V ' ^ 19) 
As can be seen from ( |3.18 ), ( 3.19| ), ( |3.16[ ), (3.15) the symmetry with respect to the change 7 to —7 inherent in 



the spin model (2.2) is destroyed in the elaborated formulas. Evidently, the derived results for the XZ Hamiltonian 
in the limit 7 = — 1 coincide with the results for the XY Hamiltonian with 7 = ±1 presented in the first part of 
this Section. 

It can be easily observed that in the limit 7=1 the described 2D Jordan- Wigner fermionization approach for 
the square-lattice s = \ XZ model yields the mean-field theory of the square-lattice Ising model. Really, Eq. 



( 3.12 ) for 7 = 1 (in this case the fermionic representation ( 3.12 ) becomes exact since the Ising interaction of z spin 
components does not involve the phase factors) after decoupling ( [3.13 ) becomes 

r 



Afe-l Ny-1 . 

# = E E ( 2J - 2 ) K+i,^+m - 2 ) + 2J± " 2 

= E E (2^ ( ms ij- + msf +lj - — m 2 ) + 2Jj_ (msjj + ms|j +1 — to 2 )) 



(3.20) 



i=0 3=0 

(we have used the relation s z = d + d— |). Obviously, the Hclmholtz free energy per site for the mean-field spin 



Hamiltonian (3.20) follows also from (3. IS) after inserting 7 = 1, i.e., / = — i In (2 cosh (2(3 (| J| + | J±\) m)) 
2 (| J | + | Jj_|) m 2 . The equation for to (3.1£) becomes simple 



to = -tanh(2/3(|J| + \ J±\)m) 



(3.21) 



and hence we arrive at familiar formulas for the internal energy per site e = — 2(|J| + |Jj_|)m 2 and the specific 



heat 



4/3 2 (|J| + |J ± |)to 



dm 
9/3 ■ 



To close this Section let us note, that a more sophisticated treatment (i.e., a more complicated decoupling of the 
l.h.s. of ( |3.13j )) of the Hamiltonian ( |3.12 ) is possible (see Ref. M for a review of analogous studies for the isotropic 
Heisenberg model). However, we found that this more complicated approximation seems to be less adequate to 
describe the properties of the spin system. 



4 Finite— size scaling analysis of the exact diagonalization data 

In this Section we present the results of the exact diagonalization computations of the ground-state energy and the 
ground-state magnetization for the s = ^ isotropic (i.e., 7 = 0) XY model on small spatially anisotropic (i.e., < 
R < 1) lattices with periodic boundary conditions. The results for eo a,ndm x2 — N 2 N2 Yli[\=i Yufi j 2 =i( s ii ji s i2 h) 
for spatially anisotropic finite square lattices of N = N x N y sites (N x = N y = \/N = L) for L — 4 and L = 6 of 
square shape are reported in Table [l]. We also calculate eo and m x2 = p Em^i Emj^^mi^) ^ or spatially 
anisotropic finite square lattices of N = 8, 18, 32 sites of diamond shape (Table ||). It is worth to remark that 
according to the adopted definition of m x2 both the results for the square and diamond shaped finite square lattices 
for R = yield the adequate value for a chain divided by the number of isolated chains into which the finite square 
lattice splits for R = (compare, for example, the results from Table [l] and Table || for m x2 (L = 4) (4 chains at 
R = 0) and m x2 (N = 8) (2 chains at R = 0) or for m x2 (L = 6) (6 chains at R = 0) and m x2 (N = 18) (3 chains at 
R = 0)). 
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Table 1: Exact diagonalization data for the spatially anisotropic s = \ XY model (7 = 0) on finite square lattices of square 
shape. 



R 


e (L = 4) 


eo(L = 6) 


m x2 {L = 4) 


m x2 (L = 6) 


0.00 


-0.353553 


-0.333333 


0.045535 


0.025463 


0.01 


-0.353586 


-0.333385 


0.047178 


0.027241 


0.02 


-0.353684 


-0.333541 


0.048908 


0.029217 


0.03 


-0.353848 


-0.333801 


0.050731 


0.031427 


0.04 


-0.354078 


-0.334164 


0.052652 


0.033911 


0.05 


-0.354375 


-0.334632 


0.054673 


0.036713 


0.10 


-0.356909 


-0.338570 


0.066284 


0.056225 


0.15 


-0.361291 


-0.345166 


0.079674 


0.078399 


0.20 


-0.367541 


-0.353772 


0.092803 


0.092460 


0.30 


-0.384698 


-0.374145 


0.111863 


0.105124 


0.40 


-0.405918 


-0.396792 


0.121758 


0.110963 


0.50 


-0.429449 


-0.420846 


0.126876 


0.114331 


0.60 


-0.454432 


-0.445916 


0.129709 


0.116406 


0.70 


-0.480432 


-0.471769 


0.131334 


0.117683 


0.80 


-0.507203 


-0.498253 


0.132243 


0.118429 


0.90 


-0.534590 


-0.525255 


0.132689 


0.118804 


1.00 


-0.562486 


-0.552694 


0.132816 


0.118911 
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Table 2: Exact diagonalization data for the spatially anisotropic s = \ XY model (7 = 0) on finite square lattices 
diamond shape. 



R 


e (iV = 8) 


e a (N = 18) 


e (N - 32) 


m x2 {N = 8) 


m x2 (N = 18) 


m x2 (N = 32) 


0.00 


-0.353553 


-0.333333 


-0.326641 


0.091069 


0.050926 


0.033695 


0.01 


-0.353617 


-0.333387 


-0.326712 


0.094262 


0.054471 


0.037441 


0.02 


-0.353809 


-0.333555 


-0.326929 


0.097420 


0.058335 


0.041851 


0.03 


-0.354127 


-0.333845 


-0.327298 


0.100523 


0.062478 


0.046982 


0.04 


-0.354570 


-0.334266 


-0.327828 


0.103553 


0.066833 


0.052763 


0.05 


-0.355134 


-0.334823 


-0.328526 


0.106492 


0.071311 


0.058967 


0.10 


-0.359648 


-0.339627 


-0.334342 


0.119459 


0.092111 


0.085177 


0.15 


-0.366554 


-0.347155 


-0.342674 


0.129211 


0.105816 


0.097591 


0.20 


-0.375259 


-0.356377 


-0.352257 


0.136131 


0.113676 


0.103995 


0.30 


-0.396250 


-0.377584 


-0.373529 


0.144342 


0.121607 


0.111009 


0.40 


-0.420165 


-0.400902 


-0.396613 


0.148511 


0.125525 


0.114920 


0.50 


-0.445817 


-0.425543 


-0.420942 


0.150783 


0.127801 


0.117335 


0.60 


-0.472609 


-0.451135 


-0.446204 


0.152087 


0.129200 


0.118862 


0.70 


-0.500218 


-0.477457 


-0.472195 


0.152846 


0.130060 


0.119815 


0.80 


-0.528446 


-0.504366 


-0.498777 


0.153274 


0.130561 


0.120375 


0.90 


-0.557169 


-0.531756 


-0.525847 


0.153484 


0.130813 


0.120657 


1.00 


-0.586302 


-0.559552 


-0.553329 


0.153543 


0.130885 


0.120738 
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From [|29[ gl[ |22| we know that the ground-state energy per site for the case R = 1, 7 = scales as 

e (L) = eo(oo) + + g- + . . . . (4.1) 

From Rcf. |3(]] we know that the finite-size corrections to the ground-state energy per site for the case R = 0, 
7 = (for both periodic and antiperiodic boundary conditions) scales as 

e (N x ) - eo(oo) ~ • (4.2) 

To extrapolate the exact diagonalization data with periodic boundary conditions for < R < 1 we use the scaling 
law 

e (L) = e (oo) + (4.3) 

which contains two unknown parameters eo(oo) and e and yields correct asymptotics in ID and 2D limits. The 
same scaling law was used in |l^] to extrapolate the exact diagonalization data for spatially anisotropic 2D s = | 
Heisenberg model. Another scaling law which also yields correct asymptotics in ID and 2D limits and therefore 
may be assumed to extrapolate the exact diagonalization data for < R < 1, 

e (L) = eo(oo) + + + . . . , (4.4) 

contains three parameters eo(oo), £2, £3 to be determined from the exact diagonalization data. Bearing in mind the 
2D-to-lD crossover we should note, that the available data for N = 16, 36 square lattices cannot be completed 
by data for N = 25 square lattice since the latter lattice for R — > transforms into chains of odd number of sites 
and the data essentially differ from those for chains of even numbers of sites. The exact diagonalization data for 
N — 64 are far beyond the available computer resources. On the other hand, the lattices of N = 8, 18, 32 sites in 
the limit R — > split into 2, 3, 4 chains of 4, 6, 8 sites and the shape effects become more pronounced. 

The simplest assumption for extrapolation of the exact diagonalization data for m x2 (L) for R = 1, 7 = (see, 
for example, [^|[ [[l]) is as follows 

m x2 (L) = m x2 (00) + — + .... (4.5) 
L 



We use the finite-size scaling (4.5) to extrapolate the numerical data for N = 16, 36 and for N ~ 8, 18, 32 lattices 
for < R < 1. The finite-size scaling analysis of the available data yields incorrect prediction m x2 (oo) < in 
quasi-lD limit starting from R in between 0.03 and 0.05 (see Section 5 below) that may serve as indication of a 
disappearance of the long-range order in ID limit. 
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From J32| we know that in the ID case R = 0, 7 = x (or y) magnetization vanishes as 

- 2 ,„, 1 



m xz (N x ) 



(4.6) 



To obtain m x \oo) for < R < 1 we may combine (|4.q) and (fhq) and assume to 1 (L) to behave as follows 



or as follows 



m- 2 (L) = m- 2 M + ^ 
L 2 



m x2 (L)=m* 2 (oo) + -^ + ^ 1 



(4.7) 



(4.8) 



The scaling which appears in (4.7) and (4.8) as i? — > (as well as in ( |Q| ) for i? = 0) when applied to accessible 
small systems works poorly probably because L is still too small. 



To summarize, using the data for N = 16, 36 and N = 8, 18, 32 clusters and scaling laws (13) and ( |4.5| ) we 



find eo(oo) and m x2 (oo) (and the order parameter m = y/m?' 1 + m,y 2 = V2m x2 ). We use these findings for the 
discussion in the next Section. Let us also remark, that the presented discourse refers to the case 7 = 0, whereas 
the case of anisotropic exchange interaction 7 7^ was not considered. 



5 Effects of the spatial and exchange interaction anisotropics 

In this Section we discuss the results for the ground-state energy, the order parameter and the specific heat obtained 
by the different methods illustrated in Sections 2, 3, 4. We also compare the approximate and exact results when 
available. 

In Fig. ^ we show the ground-state energy of the square-lattice s = | XY model obtained within various 
approaches in comparison with some exact results (the ID limit and the Ising limit). As can be seen from Eq. 
( |3.9[ ), (|3.6| ) the Jordan-Wigner fermions for the XY model (short-dashed curves) recover the ID limit (R = 0, 
R = 00). In the other limit of the Ising model (7=1) this approach fails (obviously except the ID limit R = 0, 
R = 00). The Jordan-Wigner fermions for XZ model (long-dashed curves) correctly reproduce the dependence eo 
vs. R for 7 = 1 but not the exact results for R = 0. We also show the dependence eo vs. —7 obtained within the 
latter approximation (thin long-dashed curves with cusps) to demonstrate explicitly a difference in the outcomes 
as 7 changes its sign. It should be remarked that the linear spin-wave theory predictions for the ground-state 
energy ( |2.9| ) (dotted curves) are in a satisfactory agreement with the exact results for R = (except small 7), 
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coincide with the exact result for 7 = 1 and are in a good agreement with the exact diagonalization data for 7 = 
(full squares). 



The XZ model within the fcrmionic picture exhibits in-plane magnetization introduced by decoupling (3.13) 



which is determined self-consistently by Eq. (3.19). Fig. |] shows the resulting dependence of the ground-state 



magnetization m z (long-dashed curves) on 7 for R = 1 (curves 3), R — 0.5 (curves 2), and R = (curves 1). 
The value of m z for R = 1 overestimates the exact diagonalization result at 7 = (full square). Moreover, the 
magnetization does not disappear in ID limit at 7 = that contradicts the rigorous statement and indicates the 

incorrect treatment of strong quantum fluctuations in one dimension. The exact zero-temperature relation between 

1 

2\ 4 



in-plane magnetization and exchange anisotropy parameter in ID reads m — ^ — \J~^J J IH] (solid curves 
1 in Fig. ^). The linear spin-wave theory results ( 2. 10] ) (dotted curves) are valid only for 7 > as was explained 



in Section 2. (Therefore, within the frames of the linear spin-wave theory we cannot trace the behaviour of m z 
(for example, for R = 1) when decreasing 7 to negative values.) Contrary to the fermionic description the linear 
spin-wave theory predicts smaller values of m z at 7 = in a good agreement with the exact diagonalization data 
(see the upper inset in Fig. ^) . The linear spin-wave theory value of m z for 7 = decreases as R decreases and 
finally becomes zero at a certain small but nonzero value of R. If i? — > the value of m z for 7 = tends to —00 
which is a manifestation of the inapplicability of the linear spin-wave theory in the ID limit (see the lower inset in 
Fig. H where the spin-wave theory result for R = is shown). 

Let us consider the dependence of the ground-state in-plane magnetization on 7 as it issues from the fermionic 
picture in more detail. Assume, for example, the spatially isotropic case R= 1. From exact results for this model 
we know that for 7=1 the order parameter to — ^ is directed along z axis, whereas for 7 = +0 the order parameter 
is < to < I and it remains directed along z axis. If 7 becomes negative (tending to —1) two possibilities are 
plausible: i) for infinitesimally small 7 = —0 the order parameter m turns to the x axis and its value starts to 
increase (approaching i as 7 — > —1); ii) the order parameter to remains directed along z axis and its value decreases 
until a certain characteristic value of (negative) 7 at which m abruptly turns to x axis and then its value starts 
to increase (approaching ^ as 7 — > — 1). The fermionization approach results give evidence in favour of the latter 



scenario. To illustrate this issue we show in Fig. rathe ground-state energy (given by (3.18|) as (3 — > 00) for spatially 



isotropic {R = 1) XZ model as a function of m z or m x (see the comments after Eq. (3.13)) for several values of 
7 < 0. The profiles plotted in the left panel exhibit one minimum at m z ^ for 71 < 7 < (curve 1), two minima 
at m z (the deeper one) and at m z — for 72 < 7 < 71 (curve 2), two minima at m z ^ and at m z = of 
the same depth for 7 = 72 (curve 3), two minima at m z ^ and at m z = (the deeper one) for 73 < 7 < 72 
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(curve 4), and one minimum at ra z = for 73 < 7. (The curve m z vs. 7 in Fig. ^jjumps to zero value at 7 = 73.) 
The profiles plotted in the right panel represent the dependence of the ground-state energy for the same 7s on the 
order parameter directed along x axis. From these plots we see that for 7 < the minimum in the dependence 
eo vs. m x is always deeper than the one in the dependence eo vs. m z and hence the x directed order parameter 
is favourable. However, while decreasing 7 from positive values to negative ones the system may remain in the 
metastable phase with the order parameter directed along z axis until 7 achieves 73. This is a typical scenario 
of the first-order quantum phase transition driven by the exchange interaction anisotropy 7. The corresponding 
region of metastability in the 7 — i?-plane is shown in Fig. |J. The obtained issue can be compared with the 
predictions of other approaches. The linear spin-wave theory does not work if 7 < and thus cannot be used for 
this discussion. The coupled cluster method and the correlated basis function method similarly to the 2D Jordan- 
Wigner fermionization approach suggest an analogous behaviour of the ground-state energy J^|, ^5|. Moreover, 
the correlated basis function prediction for the value of 7 = —0.36 associated with a phase transition for spatially 
isotropic square-lattice p5| is quite close to the estimate obtained by the fermionization approach 73 = —0.3528. 
An intriguing question naturally arises here: is the described scenario of a first-order quantum phase transition 
inherent in the spin model or it is an artifact of the approximate approaches? This question, apparently, cannot 
be answered by the mentioned theories and more work is necessary to draw a definite conclusion. 

Next we pass to another interesting problem, that is, the appearance/disappearance of the long-range order at 
zero temperature in quasi-lD limit. It is known that there is no long-range order for the ID (R = 0) isotropic 
(7 = 0) XY model and it does exist for the spatially isotropic 2D (R — 1) isotropic (7 = 0) XY model. If 7 = 1 
the order parameter takes its classical value m = | for all < R < 1. Let us consider the dependence m vs. R 
which arises within different approaches. The prediction for the order parameter m z vs. R as it follows from the 
Jordan- Wigner picture for the XZ model with 7 = is shown in Fig. ^| by long-dashed curve. The exact result 
m z = for R = is not reproduced. The result of the linear spin-wave theory is shown in Fig. ^| by dotted curve. 
The linear spin-wave theory does not work properly as R — > yielding m z — > — 00. To make this explicit, we insert 



(|2.7| ) into ( 2.1C ) and note that the integrand diverges at R = 0, 7 = around k x — 0. We expand the integrand 
about small k x to pick up the main contribution to the integral in the limit R, 7 — > 0. For 7 = 0, R — > we get 

m z = — =— In R + finite terms. (5.1) 
4V2tt 

For finite 7 > the order parameter m z is finite for all < R < 1; if R = the order parameter again has a 
logarithmic divergence as 7 — > 0, m z — 4v ^ 7r In 7 + finite terms; for finite R > the order parameter m z is finite 
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for all < 7 < 1. The dotted curve in Fig. ^]for very small R reproduces the logarithmic divergence for 7 = and 
the characteristic value of R at which m z becomes zero is about 7 • 10~ 4 . We also report the exact diagonalization 
data for the N — 36 lattice shown by empty squares in the main plot in Fig. || as well as the finite-size scaling 



(4.5) predictions from square and diamond shaped clusters shown by full squares and full diamonds, respectively. 



According to ( |4.5| ) m x (oo) 2 becomes negative for R in between 0.03 and 0.04 (diamond shaped clusters) or in 
between 0.04 and 0.05 (square shaped clusters) (see, for example, two lower curves in the inset in Fig. [| which 
refer to the square shaped clusters) that we may take as a criterion for disappearance of the long-range order as 
the spin system is becoming a collection of noninteracting chains. 

We may summarize our findings as follows: for the spatially anisotropic 2D s = i isotropic (i.e., 7 = 0) XY 
(XZ) model the 2D Jordan- Wigner fermionization approach predicts long-range order for any < R < 1, the 



linear spin-wave theory for R > R c with R c about 7 • 10~ 4 , the finite-size scaling data (4.5) for R > R c with R c in 
between 0.03 and 0.05. It should be remarked that the results of the performed analysis are qualitatively similar to 
what have been obtained for weakly coupled s = \ isotropic Heisenberg chains |l(| [El], [l2| (see also a recent study 
for a 2D array of interacting two-leg s = \ Heisenberg ladders ]lj|)- Note, however, that the exact diagonalization 
data prediction for R c for the XY Hamiltonian is in between 0.03 and 0.05 whereas for the Heisenberg Hamiltonian 
R c is essentially larger been in between 0.2 and 0.3 (see Fig. 2 of Ref. For both systems, i.e., the isotropic XY 

model and the isotropic Heisenberg model the question about developing of the long-range order as the interchain 
interaction starts to increase is far from being settled and more delicate approaches are necessary to resolve this 
question. 

Let us turn to the properties of the considered square-lattice s = | model at nonzero temperatures. In Fig. ^ 
we show the temperature dependence of the specific heat for the XY model ( ^3 -11 ) (short-dashed curves) and XZ 



model (which follows from (3. IS), (3.19)) (long-dashed curves) in comparison with the exact results (solid curves) 



for R = and the Onsager solutions (for 7 = 1) fl33| as well as with the quantum Monte Carlo result for 7 = Q, R = 1 



[p.9| (open squares). Eq. (3.11) yields exact result in ID limit. The noninteracting Jordan-Wigner fermions (3.7), 
( |3.6| ) cannot reproduce the logarithmic singularity inherent in the spin model for 7=1 (however, recover the exact 
result as the square-lattice Ising model becomes a set of noninteracting Ising chains). The temperature dependence 



of the specific heat for XZ model as it follows from (3.18), ( p. 19 ) has a mean-field character exhibiting an incorrect 



finite jump. A noticeable difference between the outcomes of the quantum Monte Carlo and 2D Jordan-Wigner 
fermionization approaches for 7 = 0, R = 1 is seen at low and intermediate temperatures. Recently the temperature 
dependence of the specific heat for 7 = 0, R = 1 has been examined within the Green's-function approach p6| . 
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The obtained temperature profile (Fig. 2 of that paper) although is in qualitative agreement with the quantum 
Monte Carlo data but exhibits a lower peak at a higher temperature. 

In Fig. [?] we plot the temperature dependence of the in-plane magnetization for spatially isotropic 2D (R = 1) 
XZ model derived within the 2D Jordan-Wigner fermionization scheme (long-dashed curves) and the exact result 
for 7 = 1 (solid curve). For R = 1, 7 = 1 the approximate result based on the fermionization approach (long- 
dashed curve 3 in Fig. [?]) is the mean-field counterpart of the Onsager solution (solid curve 3 in Fig. [7]). The 
approximate result for R = 1, 7 = (long-dashed curve 1 in Fig. Q) is incorrect since the Mermin- Wagner theorem 
predicts m = at any finite temperature for the isotropic XY model. For R = the fermionization approach 
yields decreasing of m z as temperature increases and vanishing of m z at a certain temperature instead of the correct 
behaviour m = for any finite temperature. 

6 Summary 

To summarize, we have extended the 2D Jordan-Wigner fermionization approach for the anisotropic square-lattice 
s = 4 XY model. In our study of the ground-state properties we have concentrated on the effects of i) exchange 
interaction anisotropy and ii) spatial anisotropy. We have analyzed 1) the quantum phase transition driven by the 
exchange interaction anisotropy and 2) appearance/disappearance of the long-range order in quasi-lD limit. We 
have discussed these issues using different approaches: 2D Jordan-Wigner fermionization, linear spin-wave theory 
and exact diagonalization data supplemented by finite-size scaling. The study of the anisotropy driven changes 
extends the the earlier papers |25| suggesting the first-order quantum phase transition scenario. The study of 
the 2D-to-lD crossover constitutes the XY counterpart of the problem intensively discussed within the context of 
the Heisenberg model (jl2| and references therein). In our study of the temperature effects we have discussed to 
what extent the Jordan-Wigner fermions can reproduce the peculiarities of the specific heat for the square-lattice 
Ising and isotropic XY models as well as the relation to the Mermin- Wagner theorem. 

We have found that the 2D Jordan-Wigner fermionization for the XY Hamiltonian works well in quasi-lD 
limit (especially as 7 — > 0) for calculating both the ground-state and thermodynamic quantities. Being applied 
to the XZ Hamiltonian the 2D Jordan-Wigner fermionization gives the ground-state results similar to the linear 
spin-wave theory. Although both approaches reproduce the existing long-range order in two dimensions (and 
therefore work well as 7 — > 1), they overestimate its effects and fail as R — > 0. The temperature effects for the XZ 
model are poorly reproduced within the 2D Jordan-Wigner fermionization scheme. 
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Figure 1: The ground-state energy of the square-lattice s = \ XY (XZ) model: exact results (solid curves), finite- 
size extrapolation (N = 16, 36) (full squares), linear spin-wave theory (dotted curves), fermionization approach for XY 
model (short-dashed curves) and XZ model (long-dashed curves). Thin long-dashed curves with cusps in the right panels 



correspond to the dependence eo vs. —7 as it follows from ( 3.10 , ( 3.19 ) in the limit —* 00. The inset demonstrates size 



scaling ( L3 ) of the finite lattice data (N = 16, 36) for R = 1 2 §.8, 0.6, 0.4, 0.2, (from bottom to top). 
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Figure 2: The zero-temperature magnetization m z vs. the exchange interaction anisotropy parameter 7 for square-lattice 
s — i XZ model (1: R = 0, 2: i? = 0.5, 3: i? = 1): exact result for R = (solid curves), finite-size extrapolation result 
(N = 16, 36) for R — 1 (full square), linear spin-wave theory results (dotted curves) and fermionization approach results 
(long-dashed curves). In the lower inset the spin-wave theory result for R = becomes negative for 7 less than 1.7 • 10~ 5 . 
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Figure 4: Region of metastability for the square-lattice s — \ XZ model as it follows from the fermionization approach. 
For all 7 < the stable phase has m z = 0, m x ^ 0. In the region A the phase with m z 7^ is still possible although not 
favourable (metastable phase), in region B the phase with m z / can not exist. The solid curve which separates the regions 
A and B represents the dependence 73 vs. 7? (see explanation in the main text). 
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Figure 5: The order parameter for the square-lattice s = | isotropic (i.e., 7 = 0) XY {XZ) model as the spatial anisotropy 
parameter R varies from to 1: exact diagonalization data for N — 36 square shaped clusters (empty squares in the main 



plot), exact diagonalization data completed by scaling (4.5) (square shaped clusters: full squares, diamond shaped clusters: 
full diamonds), linear spin-wave theory (dotted curve), fermionization approach (long-dashed curve). The linear spin-wave 



theory result for m becomes zero at 7? about 7-10 4 . The inset demonstrates the size scaling (4.5) of square shaped clusters 
for R — 0.04, 0.05, 0.1, 0.2, 0.4, 0.6, 0.8, 1 (from bottom to top; the curves for R = 0.8 and R — 1 coincide in the chosen scale). 
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Figure 6: The temperature dependence of the specific heat for the square-lattice s = 1 XY {XZ) model: exact results 
(solid curves), quantum Monte Carlo result (open squares) ju^, fermionization approach results for XY model (short-dashed 
curves) and for XZ model (long-dashed curves). The left panels refer to the case 7 = 0, the right panels refer to the case 
7 = 1. For R = (two upper panels) the solid and short-dashed curves coincide. 
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Figure 7: The temperature dependence of the order parameter m z for spatially isotropic (R = 1) s = \ XZ model (1: 
7 = 0, 2: 7 = 0.5, 3: 7 = 1): exact (solid curve) and fermionization approach (long-dashed curves) results. 
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